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CHAPTER I, 

INTRODUCTION 

In the operation of future manned space vehicles 
there is always a finite probability that an accident 
will occur which results in uncontrolled tumbling of 
a spacecraft. This uncontrolled motion creates a 
hazardous environment to the crew, which would experience 
oscillating accelerations. The structural integrity of 
the disabled vehicle may be jeopardized by prolonged 
tumbling presenting additional danger. An earth-launched 
rescue vehicle may not arrive for approximately twenty- 
four hours due to fueling, launching, and rendezvous 
operations. Hard docking by a manned rescue craft with 
the disabled vehicle is not possible because of the 
ha.zardous environment to which the rescue crew would be 
exposed and the excessive accelerations and fuel usage 
required of the rescue vehicle. Detumbling would then 
have to be accomplished by an external means using a 
remotely controlled tele-operator or by impinging fluid 
jets on the disabled craft. Therefore, it is desirable 
to develop an internal autonomous control system to 
either completely detumble the vehicle or lessen the 
tumbling motions until the rescue craft arrives. Such a 
device would become active upon loss of control. 
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Mass expulsion devices require onboard storage of 
propellant and, hence, may not be reliable for this 
application on a long term basis. Momentum exchange 
devices, such as control moment gyros, may require con- 
tinuous operation since startup of the gyros would .be 
difficult once a tumbling situation has occurred. Power 
requirements for continuous operation of this safety 
device may be prohibitive. Control moment gyros exper- 
ience a marked increase in gimbal loading in a tumbling 
vehicle, and would be subject to gimbal failures. These 
devices also have a tendency for saturation in large 
corrective maneuvers. Passive energy dissipation devices 
such as viscous ring and pendulum dampers are reliable 
and simple but have relatively low energy dissipation 
rates, and therefore, require a long time for stabiliza- 
tion. These devices seem most appropriate for vehicles 
which have a high nominal spin rate about one axis. For 
this application passive devices can effectively damp 
out relatively small transverse rates. However, an 
emergency detumbling system should be capable of stabiliz- 
ing a disabled vehicle with arbitrary tumbling motions 
where all angular velocities may be of the same order of 
magnitude . 

The development of a. movable mass. control system to 
convert the tumbling motions of a disabled vehicle into 
simple, spin is presented here. A simple spin state would 
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greatly facilitate crew evacuation and final despinning 
by an external means. The system moves a control mass, 
according to a selected control law, in the acceleration 
environment created by the tumbling motion. By moving 
the mass properly, the rotational kinetic energy of the 
system may be increased or decreased creating simple spin 
states about the minimum or maximum moment of inertia 
axis, respectively. The control system is designed for 
the latter case due to its associated stability in the 
presence of perturbing forces. 

The complete equations of motion of a rigid space- 
craft with attached control mass are formulated making no 
assumptions concerning vehicle symmetry or magnitude of 
the transverse angular rates. A control law relating 
control mass motions to vehicle motions is selected 
based on Liapunov stability theory. A method of determin- 
ing control system parameter values, based on an estimate 
of the worst case tumble state, is also presented. For 
a large space station and realistic initial conditions, 
it is shown that the movable mass control system is 
capable of decreasing the kinetic energy of the system 
to its minimum state, establishing a simple spin about 
the maximum moment of inertia axis, in several hours. 

An additional application of a movable mass control 
system as a wobble damper for an artificial-g mode, space 
station is also presented. Finally, comparisons with 
competitive passive devices for this application are made. 
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CHAPTER II 

PREVIOUS INVESTIGATIONS 

The effect of internal moving parts on the attitude 
motions, of a space vehicle has been investigated by 
Roberson (1) and Grubin (2). Both of these authors have 
developed the equations of motion of a rigid main body 
carrying an arbitrary number of moving rigid bodies. 

These studies are primarily concerned with analytically 
describing the result of internal mass motion so that 
possible destabilizing effects may be determined. 

Roberson chose the composite center of mass of the 
system as the reference point for the equations of motion. 
For this choice of reference point, the formulation leads 
to time varying moments of inertia of the main vehicle 
since the reference point is moving relative to the 
vehicle due to mass motion. Grubin circumvented this 
problem by choosing the center of mass of the vehicle as 
the, ref erence point for the equations of motion. For 
this choice, the moments of inertia of the vehicle are 
constant relative to a body fixed coordinate system. As ■ 
noted by Grubin, the formulation may produce, equations 
of motion which are no simpler than those produced by 
Roberson's method. However, the selection of the vehicle 
center of mass as the reference. point seems more natural 
since the rotational motions of the spacecraft may be 
referred directly to the vehicle center. of mass rather 
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than the moving composite center of mass of the system. 

In addition, the motion of the moving mass may be ex- 
pressed relative to the vehicle axes which is of some 
practical advantage. 

Kane and Sobala (3) investigated the possibility of 
using the effect of internal mass motion to provide atti- 
tude stabilization. The system considered was a sym- 
metric . satellite and two particles performing prescribed 
oscillations. The particles were constrained to move 
along the axis of symmetry. In a later paper Kane and 
Scher (4) discussed the possibility of moving the parts 
of a connected system relative to each other to convert 
undesirable tumbling or nutation into simple spin. The 
proposed concept exploits the following facts. When the 
system . is constrained to move as a single rigid body, the 
rotational kinetic energy is constant and bounded between 
two limits. These limits .are determined by the angular 
momentum and the maximum and minimum moments of inertia, 
of the system. If the parts- of the system are moved 
relative to each other and again constrained to form a 
rigid body, the rotational kinetic energy of the new 
configuration will, in general , be different from the 
initial value. Therefore, the control scheme presented, 
by Kane and Scher is to move the parts relative to each 
other in such a way that the -rotational kinetic energy 
changes to its maximum or minimum value. The system 
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would then be in a simple spin state about either the 
minimum or maximum moment of inertia axis. However, 

Kane and Scher do not present a general control law to 
accomplish this stabilization technique. It is suggested 
that a catalog of control laws be compiled for use with 
various tumble states. Then, once a tumble state has 
been identified, the proper control law may be selected 
for this case. 

Childs (5) developed a movable mass control system 
for an artificial-g space station. The analysis is not 
dependent on vehicle symmetry but the assumptions of a 
high nominal spin rate about one axis and small trans- 
verse rates are made. This allows Childs to linearize 
the equations of motion and formulate a relatively simple 
control law. However, the control system does not damp 
the transverse rates to zero but only to a constant value. 
This does not completely detumble the vehicle but does 
reduce the disturbance that crew members would experience. 

Beachley (6) suggested another application for 
movable masses. He has shown that the inversion of a 
spin-stabilized spacecraft may be accomplished using 
movable masses i The concept involves moving the masses 
to make the spin axis unstable by changing the moments 
of inertia of the system. This will cause the cone angle 
between the angular momentum vector and the line of 
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symmetry to increase. As the spacecraft completes the 
inversion phase the control masses are returned to their 
original position. The inverted spacecraft will again 
be spinning in a stable manner. 

Lorell and Lange (7) have proposed using two pairs 
of movable masses to provide an automatic mass trim 
system for a spinning spacecraft. For an active control 
system, the accuracy of the pointing control depends on 
knowledge of the relative location of the spin axis and 
the sensor axis. If the spin. axis shifts and the spin 
and sensor axes are not correctly aligned , a pointing 
error arises which is difficult to eliminate. Therefore, 
Lorell and Lange propose a movable mass system to trim 
the location of the spin axis to align the sensor and 
spin axes. The analyses presented is limited to the 
assumption of a symmetric vehicle and small transverse 
rates . 

These previous investigations have used movable 
masses for a variety of applications. However, most of 
these applications have been developed for vehicles 
which are symmetric or experience only small transverse 
rotation rates which permits some degree of linearization. 
Available literature indicates that a movable mass control 
system has not been developed for the general case of an 
asymmetric vehicle with arbitrarily large rotation rates 
about each axis. 
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CHAPTER III 

ANALYTICAL INVESTIGATIONS 

Development of Equations of Motion 

The complete equations of motion of a rigid space-. 

craft with an attached mass are developed in this section. 

In the following analyses v implies differentiation with 

4. 

respect to an inertial reference frame and [v] implies 
differentiation relative, to a body fixed reference frame. 

The generalized equation of motion for a system of 
connected rigid bodies is given by Grubin (2) as 

= S + ? X a (1) 

where M, H, and S denote the external moment, angular 
momentum, and first moment of mass of the system, re- 
spectively. The quantities are specified with respect 
to an arbitrary reference point moving in an arbitrary 
manner. The inertial acceleration of the reference point. 
IS a. This equation reduces to the standard form, M = H, 
when the reference point is fixed (a = 0) or is the 
center of mass of the system (? = 0). 

The system under consideration consists of a rigid 
main body and attached control mass shown in Figure 1. 

The reference point. for the equation of motion is selected 
as the center of mass of the main body. The body fixed 
coordinate system is aligned with the principal axes of 




Figure 1. Main Body and Attached Mass System Geometry 
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the vehicle, , X2 j and X^ , and have associated unit. 

/V /N >N 

vectors i, j, and k, respectively. 

With the assumption of no applied external torques , 
Equation 1 becomes 

i + $ X a = 0 (2) 

where H is the total angular momentum of the system. 

H = S, + S ( 3 ) 

b m 

The angular momentum of the main body, , is 

^ /s 

where , I2, and Ig are the principal moments of inertia 
of the . vehicle without the control mass and to^ , CO2 » a^nd 
Wg are the rotation rates about the X^ , X2 , and X^ axes, 
respectively. Differentiating this equation yields 

H, = [H, ] + 0 ) X H 
b b 

/s 

= + 0)2^2(12 - l 2 )]i 

As 

+ [l2<I)2 + - l3)]j 

+ ^ 12^3 + 0)2^0)2(12 - l 2 _)]k ( 5 ) 

The angular momentum of the control mass relative 
to the vehicle center of mass, 5 , is 

= mr X r ( 6 ) 
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It has been assumed that the control mass is essentially 
a point mass so that the angular momentum of the mass 
about its own axes may be neglected. The rate of change 
of angular momentum of the mass is then 

^ = mr X r . ( 7 ) 

m 

The term r is the acceleration of the control mass rela- 
tive to the main body center of mass. From Thomson (8), 
this acceleration may be written as 

• • M • •• 

r=(jjx (ojxr) + 03xr+2wx [r] + [rj. (8) 

The first term is the centripetal acceleration, the second 
is the tangential , acceleration, the third is the Coriolis 
acceleration, and the last term is the acceleration rela- 
tive to the body fixed reference frame. 

The first moment of mass of the system is given by 

^ ^ ^ /S /N 

S = mr = m(xi + yj + zk). (9) 


Note that the main body does not contribute to this term 
since the reference point is the center of mass of the 
vehicle . 

The inertial acceleration of the reference point is 



and, from Figure 1, may be written as 
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The term is the inertial acceleration of the composite 
center of mass of the system and is given as 



A ? 

~ M + m* 


Here.? is the resultant of the external forces acting on 
the system, M is the mass. of the -main body, and m is 
mass of the control mass. With the assumption of no 
external forces acting on the system, ?^ = 0 . Also, 
from the definition of center of mass of the system 


m 


m + M 


( 10 ) 


so that the acceleration of the reference point becomes 


m 

a = r— TT r . 


m + M 


Inserting these relations into Equation 2 yields 


H + 

“b m + M 


r X r = 0 . 


( 11 ) 


Defining the reduced mass of the system as 


A mM 
m + M 


and identifying 



( 12 ) 


as the force applied to the control mass and reacted to 
by, the vehicle, Equation 11 becomes 

= -? X .{. 


(13) 
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For prescribed motions of the control mass , this relation 
gives, in vector form, the equation of motion of the 
vehicle. It is evident that the dynamics may be con- 
sidered to be those of a rigid body being acted on by a 
reaction moment of -r x ? which is a result of mass 
motion . 

Expanding Equation 13 yields a set of three coupled, 
non-linear differential equations for the vehicle dynamics 
in terms of the vehicle angular rates (w^, (O 2 , ^.nd 

the movable mass position (x, y, z), velocity (x , y, z), 
and acceleration (x , y, z). All of these quantities are 
with respect to the body fixed principal axes (X^, X 2 , 

X^ , respectively). The equations are: 

Cl^ + y(y^ + z^)](L^ + Cl^ - ^2 ^ 

2 2 

+ yC-xyd )2 - xzcOg + (2yy + 2zz)w^ + yzCw^ - ^ 2 ) 

- 2 xyw 2 - 2 xzo )2 - xza)^W 2 + xyco^w^. + yz - zy] 

= 0 (14) 

[I 2 + y(z^ + x^)w2 + “ ^3 

2 2 

+ yC-yzwg - yx(I)^ + (2zf + 2 xx)u 2 + zx(co^ - lo^) 

- 2 yzoj 2 - 2yxo)^, - yxdi^(xi^ + yzo) 2 W^ + zx - xz] 


0 


(15) 
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Cig + iJ(x^ + + [I 2 - y(x^ - y^)](jjj^a)2 

2 2 

+ y[-zx(I)^ - zyih^ + ( 2 xx + 2yy)(i>2 + xy((jj2 - w^) 

- 2 zxo)^ - 2zy(i)2 - zyco^w^ + zx(j 02O)2 + xy - yx] 

= 0 (16) 


These equations are . valid irrespective of the physical 
mechanism whereby the control mass executes its motions. 

Similarly, the force, ?, applied to the control mass 
and reacted to by the vehicle, given by Equation 12, may 
be expanded in component form. The result is 


f^ = y[x - 2yo)2 + 2zw2 - yw^ + ZCO2 + yw^W2 

2 2 

+ ZW^Wg - x((02 + Wg)] 

f2 = y[y - 2zoj^ + - zw^ + xwg + 

2 2 

+ XW2W^ - yCwg m^)] 

fg = y [2 - 2XCO2 + 2 yOJ^ - XW2 yw^ t XWgW^ 

2 2 

+ yw2W2 - z(Wj^ + ( 1 ) 2 )]. 

The kinetic energy of the system is 
M 

T = 1 / t • t dm + 

From Figure 2 the following relations are evident. 


(17) 


(18) 


(19) 


( 20 ) 
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I = ^ 


S = S + f 

m c m 


Inserting these relations into Equation 20 yields 

M 

„ 1, -k ,1 f 1-5- X 

T = 7r(m + M)R *R + (•«■ / p*pdm + ^r • r ) 

L cc 2mm 

M 

*K- p dm + mr^ ) . 

The first term may be identified as the kinetic energy 
associated with the motion of the composite center of 
mass of the system. The second grouping of terms is the 
rotational kinetic energy about the composite center of 
mass. The final terms are zero from the definition of 
center of mass of the system. Since only attitude motions 
of the system are of interest, only the rotational kinetic 


energy will be considered further. 


1 i 

1 -^ -*■ 

T ^ = / p’pdm + T^r • r 

rot 2 J ^ ^ ■ 2m n 


( 21 ) 


Again from Figure 2 


Using these relations and Equation 10 , the rotational 
kinetic energy of the system becomes 


. 1 mM -5- 
dm + — — __ r 

2 m + M 


( 22 ) 
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The first term of this equation is the rotational kinetic 
energy of the rigid body about its own center of mass. 

The second term is the kinetic energy of the reduced mass 
about the vehicle center of mass. Therefore, the rota- 
tional kinetic energy of the system about the composite 
center of mass may be written as 

Trot = ^ * I * w + jvir • r (23) 

where I is the inertia dyadic of the rigid body. It is 
interesting to note that the rotational kinetic energy 
of the system given by Equation 23 is valid irrespective 
of the choice of a body fixed principal axis coordinate 
system. 

The expression for the angular momentum of the system 
with respect to the composite center of mass may be 
developed similarly. The result is, from Roberson (1), 

S^=I*co + yrxr. (24) 

It should be remembered that S in Equation 3 is the 
angular momentum with respect to the vehicle center of 
mass, whereas in Equation 24 is the angular momentum 
with respect, to the composite center of mass of the system. 
The equation of motion would be 



0 


(25) 
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if the composite center of mass were selected as the 
reference point for Equation 1. Since. Equation 2 5 indi- 
cates that is constant, the accuracy of a numerical 
solution of the equations of motion may be determined 
using Equation 24. 

Since the movable mass control system is to decrease 
the rotational kinetic energy of the system, the rate of 
change of this energy will be developed. Differentiating 
Equation 23 yields 

T ,=w*I*ai + r • f. (26) 

rot 

Partially expanding Equation 13 gives 

~ -V -> — -> -jt 

I*to + a)xI*ti)=-rxf 


so that 



-w • 0) X I • 


->• ->• 
0) - w 


X I. 


It may, be noted that the first term on the right hand 
side of this equation is zero. Using this relation, 
Equation 26 becomes 

T ,=-w»rxJ + r*?.. 
rot 


Cr] + 0 ) X r 


Noting that 


4 

r 
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this equation simplifies to the result 


T ^ 
rot 


[r] • ?. 


(27) 


Thus, the rate of change of rotational kinetic energy of 
the system is found to be independent of the vehicle, 
inertia properties and dependent only on the relative 
velocity of the control mass and the force applied to the 
mass. It may be noted that this result is valid irrespec- 
tive of the assumption of a principal axis coordinate 
system since Equations 13 and 23 are general. 

The equations of motion of a rigid spacecraft with 
attached control mass have been formulated in this section; 
The result is a set of three coupled, non-linear differen- 
tial equations for the rotation rates of the vehicle in 
terms of the control mass, position, velocity, and accelera- 
tion. The motion of the control mass will be specified 
by the control law which is selected. The expression for 
the angular momentum of the system about the composite 
center of mass has been presented to determine the ac- 
curacy of a numerical solution of the equations of motion. 
Since the movable mass control system is to provide energy 
dissipation, the expressions for the rotational kinetic 
energy of the system and the rate of change of this 
quantity have been developed. The derived relation for 
the rate of change of rotational kinetic energy will be 
used to develop a control, law in the next section. 
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Selection of Control Law 

Equations 14-16 determine the attitude motions of 
the spacecraft for specified motions of the control mass. 
It is the function of the control law to relate the mo- 
tions of the control mass to measurable vehicle parameters 
so that the control mass may respond to vehicle motions 
in an appropriate manner to lessen tumbling. A satis- 
factory control law should not be unnecessarily compli- 
cated and should not have excessive power or sensor re- 
quirements. It should , however, require determination 
of only measurable vehicle parameters , produce stable 
responses, and result in a final state of a simple spin 
about either the maximum or minimum moment of inertia 
axis. In the following analysis, the vehicle is assumed 
to have three distinct moments of inertia, I^ , and I^, 

and the relationship I^ > I 2 > is assigned to these 
quantities . 

By inspection of Equations 14-16, the equations of 
motion for an asymmetric vehicle with attached movable 
mass are extremely complicated due to their highly coupled 
and non-linear nature. Since the initial tumble rates 
may be large about, all three axes, the equations of motion 
cannot be simplified by linearization. However, several 
simple cases were identified and will be discussed before 
considering the general case. 
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The first special case requires that the motion of 
the mass be along a line parallel to and offset a distance 
b from the axis and passing through the X 2 axis. For 
this case, Equation 14 becomes 


[I^ + y(b^ + z^)]d)^ ■'^1-^3 ~ ^2 ^ y(b^ - z^)]u)2W2 

+ 2 yzzo)^ + ybzCcD^ - + ybz = 0. (28) 

Suppose the control law is chosen such that 

z = cw^^. (29) 

2 2 

With the assumption, >> 3yc u)^. Equation 28 becomes 


Cl, + yb^] 

“1 dbc “1 ^“3 " ‘^ 2^‘^1 


CI2 - I3 - yb^] 
ybc 


W2W3. 


(30) 


This equation indicates that for the case to^ >^29 
Equation 29 will result in damping of to zero produc- 
ing a stable spin about the. maximum moment of inertia 
axis. The control law would be easy to implement, re- 
quiring measurement of z and only. The mass would 
oscillate about its equilibrium position with decreasing 
amplitudes since would be damped. The control mass 
would return to its zero position when equals zero 
and a simple spin is reached. For the case of an arbi- 
trary tumbling spacecraft, the assumption > ^2 cannot 
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be made and Equation 29 does not provide a satisfactory 
control law for the general case. However, the result may 
be useful in designing a control system for a space sta- 
tion which has an artifical-g mode where the spacecraft 
has a large rate about one axis, say » and the control 
system is to damp out the small transverse rates and .CO 2 
Since a simple spin about the maximum moment of 
inertia axis is the minimum energy state of the system, 
it is evident that Equation 29 produces energy dissipation 
for this case. The second specialized case demonstrates 
that a movable mass control system may increase the energy 
of the system to the maximum energy state. The vehicle 
would then be in a simple spin about its minimum moment of 

inertia axis. For this case the control mass motion is to 

be along a line oriented parallel to, and offset by some 
distance a from the axis and passing through the X 2 
axis. For this . configuration Equation 16 becomes 

[Ig + y(x^ + a^)]o)2 + CI2 - + y(x^ - a^)]a)j^W2 

2 2 

+ 2yxxw2 + yax(w 2 - - yax = 0. (31) 

Suppose the control law is now chosen to be 

x = cco2. (32) 

With this choice of control law and the assumption I 2 >> 

2 2 

yc ojg, Equation 31 becomes 
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[Ig + Via ] ^ 2 2 

c5- - (I)o + 

3 yac 3. 1 23 


CI2 - + Via ] 

yac' 


Wj^W2 • 


(33) 


This equation indicates that if the product ac is chosen 
such that ac < 0 , for the special case Equation 32 

will result in damping of to zero and produce a simple 
spin about the minimum moment of inertia axis, X^. The 
properties of the control law are similar to those of 
Equation 29. 

The two preceding examples have demonstrated that 
the movable mass control system may increase or decrease 
the energy of the system and produce a simple spin about 
either the minimum or maximum moment of inertia axis 
for certain specialized cases. They also indicate that 
possibly the proper orientation for the direction of 
motion of the control mass is parallel to the desired 
final spin axis. However, since the initial tumble 
state of the vehicle is not known, the necessary assump- 
tions may not be made for the general case of tumbling. 

The control laws given by Equations 29 and 32 are there- 
fore attractive but inadequate in their simplicity. The 
formulation of a control law which will produce a simple 
spin, for an arbitrary tumble state is the subject of 
the remainder of this section. Although the movable mass 
control system could possibly force a simple spin about 
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the minimum moment of inertia axis, the control law will 
be developed to produce spin about the maximum moment . 
of inertia axis. Spin about this axis is stable in the 
presence of perturbing forces. 

The development of the control law starts with the 
theory of Liapunov stability. From LaSalle (9), the 
system of differential equation given by Equations 14-16 
and the mass . equation of motion ; produced by the control 
law will be completely stable and approach its minimum 
state if there exists a scalar function V(u) with certain 
properties. The variable u is the state vector of the 
system. The conditions which the Liapunov function must 
satisfy are 

a) V(u) > 0 for all u ^ 0 

b) V(u) < 0 for all u 

c) V(u)-»-<»as||ull->-«> 

where iiuii is the Euclidean norm of the state vector. For 
the physical system considered here, a convenient scalar 
function to use as a Liapunov function is the rotational 
kinetic energy of the system. Due to the nature of 
kinetic energy, conditions (a) and (c) will be automatic- 
ally satisfied. Thus, if a control law is- chosen such 
that condition (b) is satisfied, the system will be 
completely stable and approach its minimum state. 
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Consider the case where the control mass is re- 
stricted to move along a track parallel to the axis , 
and offset from this axis as shown in Figure 3. For this 
case, the rate of change of rotational kinetic energy 
given by Equation 27 simplifies to 


T 

rot 



(34) 


Thus, if the force applied to the control mass is 
selected as 

fg = -ycz; (35) 


Equation 34 becomes 


4 , • 2 

T . = -ycz 
rot 


(36) 


which satisfies condition (b). Using Equations 19 and 35, 
the resulting equation of motion for the mass is 


2 2 . • 

z .+ cz - (w^ + 0 ) 2)2 = ad )2 - h(I)^ - acOgW^ - (37) 


This equation has been written in a form which suggests 
that the mass dynamics are those of a second order system 
being forced by the motions of the tumbling vehicle. It 
should be noted, however, that the vehicle motions 
which produce the forcing function of Equation 37 are 
determined by differential equations which are themselves 
functions of control mass dynamics. If the effect of 
mass motion over one .cycle is small, the dynamics of the 
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vehicle are essentially those of a free uncontrolled body. 
With this assumption the forcing function of Equation 37 
may be considered to be dominated by vehicle motions 
over one cycle. This form of the mass equation of motion 
may then be used to determine the effect of the control 
law on mass dynamics. 

The control law given by Equation 35 satisfies condi- 
tion (b) and decreases the rotational kinetic energy of 
the system. Also, the forcing function vanishes when a 
final spin about the axis is established. However, 
due to the negative, decaying coefficient of the z term, 
the mass would not necessarily return to its initial 
position (z = 0). Therefore, it is desirable to modify 
the control law to ensure the return of the control mass 
to its zero position once a simple spin state has been 
reached. This would restore the full control capability 
of the control system should another tumbling situation 
arise . 

If the force applied to the control mass is modified 
such that 

^3 ■ “ViCj^z - y(c2 + (38) 

the equation of motion for the mass becomes 


= ad)o - b(I), - 


- bw2U)2 


z + c^z + C2Z 




(39) 
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The control law has been modified , to cancel the undesir- 
able negative coefficient of the z term and to provide a 
positive coefficient. With this control law, the mass 
equation of motion becomes a conventional second order 
differential equation which will provide damped motion 
of the mass and return the control mass to its zero 
position once a simple spin state is established. 

The rate of change of rotational kinetic energy 
for this selection of control law is 

T^q^ = -yc^z^ - vizz(c2 + + o)^). (40) 

Here the formulation departs from the Liapunov method 
since is not necessarily negative semi-definite. 

The f irst • term. will decrease the rotational kinetic 
energy of the system while the second term will be 
oscillatory in nature, increasing and decreasing the 
energy. If the control system , constants c^ and C 2 are 
chosen properly the secular negative semi-definite term 
will dominate over the complete mass cycle. If every 
mass cycle has a net negative value for T^^^ , the system 
will approach its minimum energy state and a simple spin 
about the maximum moment of inertia axis. From the mass 
equation. of motion, the mass would return to its zero 
position, once this state has been reached. The physical 
significance of the control law given by Equation 38 will 
be discussed in Chapter IV. 
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Selection of Control System Parameters 

For the selected control law, the rate of change of 
rotational kinetic energy is given by Equation 40 as a 
function of two control system parameters, c^ and 
The parameter values should be chosen to yield a large 
value of energy dissipation rate and, hence, a fast ap- 
proach to simple spin, subject to mass amplitude and 
power limitations. It is the purpose of this section 
to develop some guidelines for the selection of these 
parameter values . 

Before proceeding with this development some 
qualitative observations will be made concerning the 
effect of c^ and C 2 on control system performance. The 
right hand side of Equation 39 is again considered as 
the forcing function of mass dynamics so that 

z + c^z + = F (41) 

where the forcing function due to vehicle motion is 

F = ao)2 - b(I)^ - aw^w^ - bo)2W2 (4 2) 

To investigate the effect of c^ and C 2 on control 
system performance, an analogy is established between 
mass dynamics and a simple spring-mass -damper system. 

In this analogy, c^ corresponds to the damping constant 
of the damper and C 2 corresponds to the spring constant 
of the spring. 



Considering the expression for only, it may 

appear that the value of should be selected to be 
large since the negative semi-definite term is propor- 
tional to c^. However this term is also a function of 
the relative velocity of the control mass so that the 
effect of a large value of c^ on mass dynamics must be 
investigated. Referring to, the spring-mass-damper 
analogy, a large value of c^ corresponds to a strong 
damper which would limit the. velocity of the mass. 

Since the negative semi-definite term of is the 

product of c^ and the relative mass velocity squared, 
an. increase in c^ may result in a net decrease in the 
magnitude of this term. The proper value of c^ must 
be determined considering both the direct effect of 
this parameter on and the indirect , effect on , the 

mass dynamics. A similar tradeoff results when the 
parameter is considered. The second term of 
is oscillatory in nature and results in energy addition 
over part of the mass cycle. This term is necessary 
to ensure the return of the control mass to its zero 
position once a simple spin state is reached. To limit 
the magnitude of the energy addition portion of the 
mass cycle, it may appear, from Equation 40, that 
should be selected to be a small value. However, again 
considering the . spring-mass-damper analogy, a small 
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value of C2 corresponds to a weak spring which would 
allow a relatively large amplitude of mass oscillation. 
These general considerations suggest that an optimum 
set of control system parameters may exist which maximize 
the energy dissipate rate while limiting the mass 
amplitude to a specified value. The remainder of this 
section presents some quantitative guidelines to aid in 
the selection of these parameter values. 

An analytic solution of the mass equation of motion 
would provide control mass position and velocity his- 
tories as functions of the parameters c^ and C2» Using 
these solutions, the rate of change of rotational kinetic 
energy could be maximized subject to the condition of a 
selected maximum mass amplitude. This optimization could 
be performed using a Lagrange multiplier technique. 
Clearly, this is not possible since the forcing func- 
tion, F, contains angular rates and accelerations of 
the vehicle. The solution of the mass equation of 
motion would therefore require the simultaneous solu- 
tion of the . vehicle equations of motion. This system 
of differential equations does not readily lend itself 
to an analytic solution due to the coupling and non- 
linearities of these equations. On the other hand, a 
completely numerical solution of the system of equations 
provides no analytic information concerning the effect 
of the control system parameters. Therefore, a partially 
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analytical , partially numerical method was developed to 
provide some insight into the effect of the parameters 
on the performance of the control system. 

The approach adopted is based on the assumption 
that the net change in the rotational kinetic energy 
over one mass cycle is small. With this assumption, 
the dynamics of the system are approximately those of 
a free, uncontrolled vehicle over one cycle. For a 
selected tumbling state, the free dynamics may be solved 
numerically. From results of the uncontrolled case, the 
forcing function of control mass dynamics, F, may be 
constructed using Equation 42. Although F will change 
as the control system reduces the tumbling, the ini- 
tial tumble state will provide .a design point which 
may be used to size control system parameters. 

The nature of the forcing function may be investi- 
gated by considering the analytic solution of the free 
vehicle dynamics. From Synge and Griffith (10), the 
angular rates of a free asymmetric body may be expressed 
in terms of Jacobian elliptic functions. The result is 

= Y cn [p(t - tQ>] 

032 = 3 sn 


o)g = a dn [p(t.- tp)] 
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2 

for the condition H > 2I2T which corresponds to motxon 

about the maximum moment of inertia axis, . For the 
2 

condition H < 2I2T, which corresponds to motion about 
the minimum moment of inertia axis, X^ , the solution is 

0)^ = y’ dn [p'(t - tp)] 

( 1)2 = 3* sn [p'(t - t^)] (44) 

cOg = ot ' cn [p’(t - tp)]. 

The amplitudes, ot , 6, and y, and the precession fre- 

2 

quency , p, for the case H > 2I2T and the corresponding 

2 

quantities for the case of H < 2I2T may be found in 
Synge and Griffith (10). The expressions for these 
quantities are not of direct concern here since it is only 
the form of the forcing function which is being investi- 
gated. Using these elliptic solutions, the forcing 
function will be in one. of the following forms: 

F = a(p3 - aY)cnCp(t - tQ)]dn[p(t - t^)] 

(45) 

+ b(pY - a3)snCp(t - tQ)]dn[p(t - t^)] 
for > 2 I 2 T, and 

F = a(p'3' - cx'Y’)cnCp*(t - tg)]dn[p'(t - tg)] 

(46) 

+ b(p'Y'k^ - a'3')sn[p*(t - tg)]cn[p'(t - tg)] 

for < 2 I 2 T. 
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The parameter k is the modulus of the elliptic functions. 

The functions sn u and cn u oscillate between 1 and -1 

at the precession frequency p. The function dn u 

2 is 

oscillates between 1 and (1 - k ) at twice the pre- 
cession frequency. 

From these considerations, it is evident that F 
will be oscillatory in form. Therefore, using the tab- 
ulated form of F obtained from the uncontrolled case 
results, a Fourier series may be fit to the data points. 
The forcing function expressed as a Fourier series 
will provide an analytical expression which can be used 
to solve the mass equation of motion. This solution 
will provide information concerning the relationship 
between control parameters and system performance. 

The forcing function of the mass equation of motion 
expressed as a Fourier series is 


F = A, 


00 

+ z 

n=l 


f. 2mrt . o 

(A cos + B 

n T n 


2nTTt,. 
sin ) 


where t is the period of the function F. Since the 
forcing function, expressed as either Equation 45 or 46, 
contains no secular terms, Ag = 0. Also, since F is a 
smooth function, the Fourier series may be terminated 
with a finite number term with sufficient accuracy. 

The Fourier expansion of the forcing function then be- 


comes 
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r- V 2nirt . „ , 2niTt. 

F = S (A cos + B sin — ) . 

n=l ^ T . n T 

The Fourier coefficients, A and B , are determined 
numerically using a standard routine which generates 
Fourier coefficients for a tabulated function. 

The Fourier expansion may be rewritten as 

m ^ 

F = HD sin (s t + tan" 

T n n B 

n=l n 


where 


and 


. = r 

n n 


2 2 
+ B 

n 


2niT 


n 


The equation of motion of the mass then becomes 
z + c-z + c„z = HD sin (s t + tan” 5^). 

S- ^ 'I n D 

n=l n 

The particular solution of this equation may be readily 
obtained as 


m 

J = H E sin d> 
P n=l ^ 


(47) 


where 


and 


E = 


n 


y(C2 - 


d) = s t + tan 
^n n 


, A 
-1 n 

F“ 

n 


- tan 


, c, s 
-1 In 


^2 ” ®n 


T- 


Differentiation yields the control mass velocity as 
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z 

P 


m 

Z E s cos (j) . 
, n n 


These solutions for control mass position and velocity 
histories may be used to determine the rate of change 
of rotational kinetic energy. The result is 


T 

rot 


m 

-U ^ 

j = l 


m 

Z [E .E, s, {c, s . 
k=i 3 ^ ^ 


cos 


cos (j>j^ 


» 2 )sin (J> . cos • 


(48) 


Noting that 

cos (() . cos (|), = i[cOS(4l. - (J), ) + cos((|). + <j), )] 
D K z 3 Jc , 3 K 

sin <{) . cos = j[sin(4). + (|)^) + sin(ij). - 


the secular part of Equation 48 may be identified as 
the terms for which j = k. Thus, the secular part of 
the rate of change of the rotational kinetic energy may 
be written as , 


T 

sec 


u 2 2 

Z c,E ^s 
2 T 1 n n 
n=l 


Rewriting this equation to show the dependence of the 
control system parameters yields 


T 

sec 


9 9 9 

.U? " Bn 

^n=l [(c» - s + (c, s )^] 

2 n In 


(49) 
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which may be used to determine the effect of control 
system variables on the secular part of the energy 
dissipation rate. This equation indicates that an in- 
crease in. the weight of the control mass, and hence an 
increase in the reduced, mass , y, will linearly effect 

Increasing the Fourier coefficients A and B 
s 0 o ■ n n 

will increase T quadradically . This may be accomp- 
lished by increasing the amplitude of the forcing func- 
tion which corresponds to increasing the mass offset 
distances a and b. Thus, the control mass' track should 
be placed at the maximum allowable distance from the 
center of mass of the vehicle. 

The effect of c^ and C 2 is more difficult to obtain 

2 

but it is evident that c« should be of the order of s 

i n 

With this assumption Equation 49 becomes 


« 

T 

sec 



m 


I (A 
n=l 


n 


+ 



) 


and Equation 47 becomes 


z 

P 


- m A ^ + 

i- E 

^1 n=l 


B 


n 


sin (j) 


n 


(50) 


(51) 


Equation 50 indicates that c^^ should be selected to be 
a small value. However, from Equation 51 decreasing c^ 
results in increasing the maximum amplitude of the mass 
motion. These ..observations agree with the conclusions 
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drawn when the spring -mass-damper analogy was discussed. 
It can therefore be concluded that c^ should be the 
smallest value which limits control mass amplitude to 
its maximum allowable value. 

Equations 47 and 49 may be used to generate a 
nomograph for selection of control system parameters, 
c^ and C 2 * Once the control mass weight, track position, 

and an estimate of initial tumble state are obtained, 

• 

T may be calculated for various values of c, and c_. 
Using Equation 47 the corresponding maximum mass ampli- 
tudes may be determined. Figure 4 shows a typical nomo- 
graph which will result. This particular nomograph was 
generated for an example case which will be discussed 

in Chapter IV. The figure shows a family of curves for 

« 

Tsec ^ function of the control system parameters c^ 
and C 2 « Also shown in the nomograph is another family 
of curves which connect sets of c^ and C 2 values which 
result in a particular mass amplitude. As shown in 
Figure 4, each selected maximum mass amplitude has an 
associated optimum set of c^ and C 2 values which corres- 
ponds to a maximum energy dissipation rate. Thus, once 
the maximum allowable mass amplitude has been determined 
the proper values of c^ and C 2 may be selected from 
the generated nomograph. 
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The control system parameter selection procedure 
is simplified if the main vehicle is axially symmetric. 
For. a symmetric vehicle the modulus of the elliptic 
functions, k, is zero and the Jacobian elliptic func- 
tions may be replaced with trigonometric functions . 

For this case 


sn Dp(t - t^)] = sin p(t - t^) 
cn Cp(t - t^)] = cos.pCt - t^) 
dn [p(t - t^)] = 1. 

The forcing function of mass dynamics then becomes 

F = a(p3, - ay) cos p(t - tg) + b(py - aB) 
sin p(t - tg) 

or 

F = D sinCpt + 6) 

where 

D = \J a.^ {-pQ - ay)^ + b^(py - a3)^ 


6 


tan 


-1 a(pg - gy). 
bCpy - ag) 


ptQ- 


The equation of motion of the mass may then be written as 
z + c^z + - D sin (pt .+ 6). 
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The particular solution of the mass equation of motion 
is then obtained as 


z 

P 


D 


+ (c^pV 


sin (pt + 6 


- tan 


-1 


c,p 


1 ) 
7^- 


Co - P 


(52) 


By noting the similarity between Equations 47 and 52, 
the secular part of the rate of change of rotational 
kinetic energy may be determined immediately from 
Equation 49 by setting m = 1 and s^ =,p» Performing 
this yields 


T 

sec 


y 

2 


[ (C2 


c^D 


2 2 
P 


2 '2 X ( 7^* 

- p ) + . (c^p) ] 


(53) 


The maximum mass amplitude may be readily obtained from 
Equation 52 as 


z 

max 


D 

y<C2 - P^)^ + (c^p)^ 


(54) 


The values of c^ and C 2 which yield the maximum energy 
dissipation rate with the side condition of a selected 
maximum mass amplitude may be determined using a Lagrange 
multiplier technique. 
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Before proceeding with this formulation. Equation 41 
will be written in standard second, order form as 

z + 25a)4 + u)^z = F 
n n 


where ? is the damping ratio and is the natural 
frequency of the system. Comparing this standard form 
with Equation 41 the following relations are evident. 


Cl = 


= 2 = % 


Applying these relations to Equation 53 and rearranging 
yields 


yD‘ 


sec 


2p . .w 


CO 


n^4 


CO, 


ns 2 


+ 2(2^^ - 1)(^)^ + 1 
P P 


A 0)^ 

Defining the parameter, v = — , this equation becomes 


« 

T 

sec. 


yP^ Cv 

+ 2(2C^ - l)v^ + 1 


Similarly, the maximum mass amplitude becomes 

^ - - _E 

max ^ /— t 5 7 ) ’ 

P yv + 2(2C^ - l)v + 1 

The maximum value of T corresponding to a selected 
maximum mass amplitude is determined by constructing the 


function 
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Taking the partial derivatives of G. with respect to 
V, and X, setting these expressions equal to zero and 
solving the resulting system of' equations yields the 
optimum values of ? and v. Performing these operations 
yields three values for v. 


V = 0 , V = ±1 

Obviously, v = 0 corresponds to a minimum, since T = 0 
forv = 0, and v = -1 corresponds . to a maximum which results 
in an increase in rotational kinetic energy. The desired 
value as V = 1 which corresponds to 


or 
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Thus , the natural frequency of the control system should 
be selected to equal the driving frequency of the forcing 
function for maximum energy dissipation. For v = 1, the 
maximum mass amplitude is 

D 

z = — K— 
max o 2 ^ 

2p C 
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and, hence. 


C = 


2p^z 


max 


Thus 


parameter should be chosen such that 


D 

pz 


(56) 


max 


For these values of c^ and c.^ 


T. ^ 

sec 


yDpz 


max 


From this equation, the energy dissipation rate may be 
increased by increasing the control mass weight, moving 
the mass track away from the center of mass which in- 
creases D, or by allowing the mass. a larger amplitude. 
These observations agree with the results obtained for 
the general case discussed previously. As can be seen 
the selection of control system parameters c^ and c,^ is 
simplified for a symmetric vehicle. Once the initial 
tumble dynamics have been analyzed, the parameter values 
may be calculated directly using Equations 55 and 56 in- 
stead of generating a nomograph as was done for the 
general case. 

Although the methods presented for the selection of 
control system parameters may be. in error due to the 
assumptions and approximations made, the procedures out- 
lined will provide "ballpark" estimates of c^ and 
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values. Final selection of and values must be based 
on the actual dynamics of the system obtained by solving 
the complete vehicle and mass equations of motion. 


Sensor and Power Requirements 

The sensor and power requirements of the movable mass 
control system with the selected control law will be dis- 
cussed briefly in this section. 

For the control law expressed as Equation 38 it is 
evident that the following quantities must be sensed: 
fg, the force applied to the control mass; z, the control 

mass velocity relative to the vehicle; z, the control 

2 2 

mass position; and the combination (to^ + 0 ) 2 ). The 
quantity f^ may be determined using a linear accelerometer 
mounted on the control mass. If the accelerometer is 

/s 

mounted to detect the k-component of mass acceleration, 

the sensed acceleration will.be proportional to f^. 

The mass position and velocity may be easily sensed using 

any of a number of simple devices. The combination 
2 2 

(oi)^ + ii >2 ^ niay be determined with a linear . accelerometer 

mounted on the axis. The acceleration of a fixed 
point, P, located by a vector 3 from the vehicle center 
of mass is 


3p = (0 X (w X 3) + to X 3 


to X 
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If 2 


dk the k-component of ap 
(ap)2 = 


is 

X 2v 

+ 0)2)* 


Thus, an accelerometer placed at this point and oriented 

to. sense the k-component of the acceleration at that 

2 2 

point iWill. sense a quantity proportional to , + 0 ) 2 )* 

The sensor requirements for the implementation of the 
control law appear to be modest. 

The force requirements of the control system may be 
determined by Equation 19 since this is the force applied 
by the control system to the control mass. The instant- 
aneous power required by the control system is, by defini- 
tion, the force applied to the mass times the relative 
velocity. This definition is equivalent to Equation 27 
which gives the rate of change of rotational kinetic 
energy of the system. As noted in the second section 
of this chapter T^^^ will be oscillatory, the negative 
portions of T^q^ resulting in energy dissipation. Thus, 
the power input to the mass corresponds to the positive 
portions of In fact, if a power generation system 

would be implemented with the control system, energy 
could be stored during the energy dissipation portions 
of the mass cycle and be used to power the positive T^q^ 
portions of the cycle. This possibility is not considered 
further since the primary function of the control system 
is to aid in crew rescue and not generate power. For 
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this case, the energy input required is just the integral 
of the positive portions of the curve. The force, 

power, and energy requirements for two example cases 
will be discussed in the next section. 



CHAPTER IV 


RESULTS 

To demonstrate the feasibility of a movable mass 
control system using the control law given by Equation 3 
the space station shown in Figure 5 was chosen as an 
example vehicle. The geometric axes shown in Figure 5 
are assumed to be the principal axes for demonstration 
purposes. The properties of this Modular Space Station 
(MSS) are given below. 

= 5.15 X 10^ kg-m^(3.80 x 10^ slug-ft^) 

1 2 = 6.28 X 10^ kg-m^(4.63 x 10^ slug-ft^) 

1 3 = 6.74 X 10® kg-m^(4.97 x 10® slug-ft^) 

M = 9.98 X 10^^ kg(2.20 x 10® Ibm) 

The space station shown in Figure 5 was selected 
since it is a relatively large vehicle and does not 
have an artificial-g mode. Thus, a passive damping 
system such as a viscous ring or pendulum damper would 
not be practical for this application. 

The tumble state used for this example was obtained 
from Kaplan (11). The initial rates are 

003 ^( 0 ) = -2.86 X 10”*^ radians/sec 




Figure 5. Modular Space Station Configuration 
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11 ) 2 ( 0 ) = -0.199 radians/sec 
d)„(0) = 0.103 radians/sec. 

O 

These values are based on a realistic worst case analysis 
of the tumble state produced by a collision of a shuttle 
orbiter and the MSS. Since the resultant tumble state 
is dependent on the amount of kinetic energy trans- 
ferred by the orbiter to the MSS, the initial tumble 
conditions are given by Kaplan for 10, 50, and lOO per- 
cent kinetic energy transfer. Although the 100 percent 
case corresponds to a perfectly elastic collision which 
would not be the case, it does produce the highest tumble 
rates and was, therefore, selected to test the movable 
mass control system. 

For an uncontrolled vehicle, the selected initial 
conditions result in the tumble state shown in Figure 6. 
As can be seen, this tumble state is fairly severe and 
by no means lends itself to the assumption of small 
transverse rates since the amplitudes of the angular 
velocity oscillations are all of the same order of 
magnitude . 

The control mass track is placed at the farthest 
allowable point from the vehicle center of mass and 
oriented parallel to the maximum moment of inertia axis 
as shown in Figure ,5. For this positioning, the mass 
track offset distances are 



51 



( oas/supxp^a) A:j.TOOT;aA aPinSuy 





52 


a = 13.7 meters (45 feet) 
b = 5.5 meters (18 feet). 

For these offset distances and the uncontrolled 
dynamics shown in Figure 6 the forcing function of the 
mass equation of motion, given by Equation 42, may be 
constructed. The result is shown in Figure 7. It is 
interesting to note that, comparing Equations 19 and 42, 

F is the negative of the force on a unit mass attached 
to the tumbling vehicle, with coordinates (a, b, 0). 

Thus, Figure 7 shows a typical acceleration profile that 
a crewman would experience inside the tumbling vehicle. 

As shown in Figure 7 the acceleration would be constantly 
changing magnitude and direction and would create an 
extremely hazardous environment for the crew. A simple 
spin, on the other hand, produces constant forces with 
constant directions which would greatly aid the crew in 
evacuation operations. 

Using the forcing function shown in Figure 7 , the 
coefficients of the Fourier series representing the 
forcing function were determined. Then, using Equations 47 
and 49, the secular part of the rate of change of the 
rotational kinetic energy and the maximum mass amplitude 
were determined for various values of c^ and c^. The 
resulting nomograph was discussed in Chapter III and is 
shown in Figure 4. A control mass of 99.8 kg (200 Ibm) 
which , corresponds to 0.1 percent of the vehicle weight 




on 
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was used to generate the nomograph. Different values 
of control mass weight will shift the curves up or down 
but will not affect the shape of the curves. The dotted 
lines on the nomograph gives the secular rate of change 
of rotational kinetic energy as a function of control 
system parameters c^ and c.,^. From Figure 4, smaller 
values of c^ give larger values of Also shown 

on the nomograph by solid lines are points which result 
in the same maximum mass amplitude. For each selected 
maximum allowable mass . amplitude there appears to be an 
optimum set of c^ and C 2 which results in the optimum 
value of T . Thus, once the maximum mass amplitude 

S 60 

has been determined, the proper control system parameters 
may. be chosen using the nomograph. This nomograph is 
dependent on the initial tumble state selected as the 
design point. 

Although the nomograph shows relatively large energy 
dissipation rates for C 2 = 0, the assumptions under which 
the nomograph was generated must be kept in mind. It 
was assumed that the forcing function would be purely 
oscillatory so that it could be approximated by a Fourier 
series. For the actual case, however, the forcing func- 
tion will not be purely oscillatory, in general, since 
the system is being damped. A value of C 2 = 0 corresponds 
to having no spring in the spring-mass -damper analogy dis- 
cussed previously. Therefore, the non-oscillatory nature 
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of the forcing function will cause the control mass 
oscillations to migrate away from its initial position 
for C 2 = 0. Thus, a non-zero value of C 2 is required to 
prevent this and insure the return of the control mass 
to its initial position. 

Once the control system parameters have been 
selected, the mass and vehicle equations of motion may 
be solved. This was done using a fourth-order Runge- 
Kutta algorithm to numerically solve the system of dif- 
ferential equations with an IBM 370-165 computer. The 
accuracy of the algorithm was checked using the fact 
that the angular momentum of the system about the com- 
posite center of mass is constant in the absence of 
external moments. 

Before stating the results, the physical signifi- 
cance of the control law given by Equation 38 will be 
discussed briefly. Figures 8a and 8b show typical cycles 
of control force acting on the control mass, f^j a.nd 
velocity relative to the body fixed axes, z, respectively. 
Comparing Figures 8a and 8b it is evident that the 
control law causes the control force to be generally 
opposite to the relative velocity. Figure 8c shows 
the rate of change of rotational kinetic energy of the 
system which is the product of total force and relative 
velocity. It can be seen that T, 


is generally negative 



meters/sec) z (meters/sec) 





which corresponds to energy dissipation. The positive 
portions of correspond to energy addition. These 

portions are due to the second term of Equation 40 which 
was selected to insure that the control mass oscillates 
about its zero position and returns there after simple 
spin is established. In general, energy addition is 
not desired but from Figure 8c it can be seen that 
over a complete cycle there is a large energy dissipa- 
tion and only a relatively small energy addition over 
each cycle. Therefore, it was determined that this 
energy addition was acceptable to insure return of the 
control mass to its zero position. The situation is 
further clarified by Figure 9 which shows typical cycles 
of control mass position over the same time period as 
the quantities shown in Figure 8. Superimposed on the 
mass cycle. is a schematic of the control mass with the 
directions of the mass velocity and the force acting on 
the mass shown. It is evident that the force given by 
the control law is generally a retarding force. In 
the first half cycle the mass is moving "up" and the ap- 
plied control force is directed "down" which produces 
energy dissipation. Eventually this force will overcome 
the velocity and the force and velocity will be in the 
same direction resulting in energy addition as shown in 
Figure 8c. This situation occurs only over a small part 
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Figure 9. Typical Cycle of Control Mass Position 




59 


of the mass cycle, however. After this portion of the 
half cycle, the force and velocity will again be directed 
oppositely producing more energy dissipation. A similar 
process occurs during the second half cycle and produces 
a net decrease in the rotational kinetic energy of the 
system over the complete mass cycle. A net decrease in 
the energy over every mass cycle allows the system to 
approach its minimum energy state producing a simple 
spin about the maximum moment of inertia axis. 

If a maximum mass .amplitude of three meters is 
selected, the proper values of c^ and are, from 
Figure 4 

c^ = 3.2 radians/sec 

2 

C 2 = 0.02 radians/sec • 

For a control mass weight of 998 kg (2,200 Ibm) , the 
stated initial conditions , and the selected parameter 
values, the resultant angular velocity histories are 
shown in Figure 10. Since it is not the oscillations 
themselves but the decay of the oscillation amplitudes 
which is important, only the envelopes formed by the 
oscillations are shown, in Figure 10. The control system 
effectively collapses the o)^ and envelopes to zero. 

As this is being done the mean value of the oscilla- 
tion increases to its steady spin value as it must do 
to keep the angular momentum of the system constant. 
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Thus , the control system eliminates the transverse angular 
rates and produces a simple spin about the maximum moment 
of inertia axis. 

Figure 11 shows the envelope formed by the control 
mass oscillations. Figure 11 indicates that the maximum 
mass amplitude exceeds the predicted value of three 
meters slightly. This may be attributed to the homo- 
geneous solution of the mass equation of motion which 
was neglect. in determining the nomograph shown in 
Figure 4. The results indicate that this overshoot 
occurs only during the first mass cycle and, since the 
overshoot is small, the neglect of the homogeneous solu- 
tion is justified. It may also be noted that during the 
first part of the control period the mass oscillations 
are not symmetric about the zero position. This may 
also be attributed to the homogeneous solution since 
the oscillations become symmetric as the homogeneous 
solution decays. The homogeneous solution corresponds 
to the transient response of the control system. 

Figure 12 shows the time history of the rotational 

! 

kinetic energy of the system. The figure indicates that 
the rotational kinetic energy is oscillatory due to the 
slight energy addition during each cycle but due to the 
dominant energy dissipation port ion . of the mass cycle 
the system approaches its minimum energy state. 
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Figures 10, 11, and 12 indicate that for the stated 
initial .conditions , a maximum mass amplitude of approxi- 
mately three meters, and a control mass weight of 998 kg 
(2,200 Ibm) , a movable mass control system using the con- 
trol law given by Equation 38 is capable of converting the 
tumbling motions of the MSS into simple spin within two 
hours. To investigate the effect of various parameters 
on control system performance a time constant will be 
defined. The time constant, t , is defined as the time, 
required for the control system to collapse the to^ or 0)2 
oscillation envelope to a value of 1/e times its initial 
value. Since the general shape of the performance curves 
corresponding to other parameter values will, be similar 
to those shown in Figures 10-12 this time constant will 
be used to compared control system performance for various 
cases . 

The effect of control mass weight on system per- 
formance is shown in Figure 13. The figure shows that 
an increase in control mass weight has a marked effect 
on the system time constant. The shape of the curve is 
as expected since an extremely small mass produces very 
little enerjgy dissipation while an extremely large mass 
produces a large energy dissipation rate and, hence, a 
small time constant. However, the peak power and force 
required also increases with increased control mass weight 




Figure 13. 


as shown in Figures 14 and 15. Thus, the control mass 
weight should be selected as large as possible consistant 
with :the imposed weight, power, and force limitations. 
Surprisingly, the total energy required to operate the 
control system does not vary appreciably with control 
mass. weight. The energy requirements for the example 
considered here are shown, in Table I for various con- 
trol mass weights. Thus, for this case the control 
system seems to damp the system quickly, requiring 
high values of power and force, or damp the system 
slowly, requiring lower values of power and force, such 
that the total energy required is approximately constant. 

All of the cases discussed so far have been for 

control system parameters of = 3.2 sec”^ and C 2 = 

-2 

0.02 sec and a corresponding maximum mass amplitude of 
3.7 meters (12 feet). To demonstrate the effect of 
increased mass amplitude on control system performance 
two other sets of c^ and c^ values were selected from 
Figure 4. The result is shown in Figure 16 which gives 
the variation of the defined time constant, t , with 
maximum mass amplitude for a control mass of 499.0 kg 
(1,100 Ibm). Table II gives the values of c^ and 
that were used along with the resulting maximum amplitude 
during the first and second mass cycles and the predicted 
values given by Figure 4. As can be seen from the table. 
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Table I 

Variation of Required Energy with 
Control Mass Weight 


Control Mass Weight 

Required Energy 

199.6 

kg 

(440 Ibm) 

63.9 watt-sec 

499 . 0 

kg 

(1,100 Ibm) 

68.7 watt-sec 

725.8 

kg 

(1,600 Ibm) 

71.6 watt-sec 

998 . 0 

kg 

(2,200 Ibm) 

74.6 watt-sec 








Table II 
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the , predicted maximum amplitude falls between the actual 
maximum amplitudes during the first and second mass 
cycles. It may also be noted that the larger the mass 
amplitude the larger the overshoot relative to the pre- 
dicted value indicating the increased effect of the homo- 
geneous solution. 

An additional application of a movable mass control 
system as a wobble damper was investigated. The example 
vehicle used is the NASA 21 Man Space Station which is 
a hexagonal shaped, spin stabilized, artificial gravity 
space vehicle shown in Figure 17. The properties of 
this vehicle are given below 

= I2 = 1.42 X 10^ kg-meter^ (10.5 x 10® slug-ft^) 

Ig = 2.03 X 10^ kg-meter^ (15.0 x ,10® slug-ft^) 

M = 6.21 X 10^ kg (1.37 x 10® Ibm) 

Nominal spin rate = 0.314 radians/sec (3 rpm) 

This vehicle was the basis of a study by TRW Space 
Technology Laboratories (12) to investigate the feasi- 
bility of several passive damping devices for use as 
wobble dampers. The devices studied were. a controlled 
damping pendulum, a viscous ring damper, a pendulum 
damper, and a Naval Ordinance Test Station (NOTS) mercury 
ring damper. Selection of this vehicle as an example 
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vehicle will permit comparison of the movable mass system 
performance with the competitive devices identified by 
TRW. 

The best performing passive damper identified by 
TRW is the controlled damping pendulum using three- pend- 
ulum units, each with a tip mass of 272 kg (600 Ibm) . 

To match control system properties as closely as possi- 
ble, the control mass weight for the movable mass system 
was selected as 816 kg (l,800~lbm). The wobble state 
considered by TRW was a wobble angle of five degrees 
which corresponds to the following initial conditions. 

w^(0) = 0.0391 radians/sec 

tO2(0) = 0 

WgCO) = 0.314 radians/sec 

The location of the control mass track is selected 
as shown in Figure 17. The values of track offset 
distances are 

a = 19.8 meters (65 feet) 

b = 0 

Based on these track offset distances and the stated, 
initial conditions the. control system parameters may be. 
determined by the procedure given in the third section 
of Chapter III. From an analysis of the uncontrolled 
vehicle dynamics, the amplitude and frequency of the 
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forcing function of the mass- equation of motion are 
determined as 

D = 0.136 radians/sec^ 
p = 0.136 radians/sec 

Then using Equations 55 and 56 the control system con- 
stants may be determined. 

2 

c.^ = 0.0185 radians/sec 

c, = ■ radians/sec 

1 z 

max 

The performance of the control system was determined 
for maximum mass amplitudes of three-six meters, which 
were selected as typical values. Decay of the angle 
between the angular momentum vector of the system and 
the axis of symmetry of the vehicle is presented in 
Figure 18 for these cases. The wobble angle, 0, is 
defined as 

The values of c^ corresponding to the selected maximum 
mass amplitudes and the resulting actual maximum ampli- 
tudes are summarized in Table III. It is interesting 
to note that the transient dynamics of the control 
system causes the mass to undershoot the predicted values 
of maximum amplitudes for these cases. 




Figure 18. Performance of Control System, as. a Wobble Damper 
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Table III 

Predicted and Actual Maximum Mass Amplitudes 
for NASA 21 Man Space Station Simulation 


c 


1 


(radians/sec ) 


z 

max 

(meters ) 

Predicted Actual 


0.333 

3.0 

3.0 

0.250 

4.0 

3.9 

0.200 

5.0 

4.7 


0.167 


6.0 


5.4 
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The time constant, is defined for a symmetric 

vehicle as the time required to reduce the wobble angle, 

0 , to a value of 1/e times its initial value. The re- 
sultant time constants, peak power, force, and energy 
requirements are summarized in Table IV for the cases 
considered. Also shown in the table are the time con- 
stants obtained by TRW for the passive devices they 
considered. With the exception of the case of = 0.333 
(z =3 meters), the movable mass control system pro- 
duces faster stabilization than the passive damping 
devices identified by TRW. As shown in Table IV, the 
power, force, and energy requirements of the control 
system are modest. Thus, the feasibility of a movable 
mass control system as a wobble damper for a large 
artificial-g mode space station has been demonstrated. 



Wobble Damping Performance of Movable Mass and Competitive Systems 


79 


bX) 

U 

0) ^ 
c o 
w <u 
to 
T) I 
<U -M 
U 4-> 
*H ttJ 

' 

<u 


I — I 


LO 

J- 


co 

cn 


o 

to 

CO 
1 — 1 


<U 

•M 

<u ^ 
to 

c 
o 

-p 
5 


<u 
o 

O (U 
P-. c 


03 

<U 

Oh 


CD CM rH rH 

• • • • 

CM r-' csi r^ 

iH iH CM CM 

f-H I — I fH iH 


T3 

0) 

•H 

O* 


<D 

0 i to 
■p 

O 

00 

CM 

LO 

U -P 

• 

• 

• 

« 

O (tJ 

<N 



iH 

S 5 

O 

Cu 




iH 


(tJ 

<u 

0^ 


o 

H 


to 

tu 


4-» 






-X 



a?* 

13 


1 — 1 


t — 1 

CD 

00 


CO 

O 

d 


• 

• 

• 

• 

•• 

• 

• 

• 

•H 


LO 

CO 

oo 

CM 

CO 


LO 

o 










o 










iH 


e 








ti 


tu 





bO 





p* 





a 





to 





•H 





>1 





a 


fn 



CO 





e 


<U 


tu 






ttJ 


a 



CO 


o 

o 


Q 

bO 

e 

>, 

•H 

to 

00 

LO 

o 

CD 


C 

Id 

u 

> 

(tJ 

00 

CM 

CM 

iH 

T3 6 

•H 

Q 

d 

tu 

s 

* 

• 

• 

« 

<u 

Pi 


o 

Q 


o 

o 

o 

O 

iH r— 1 

u 

S 

Pt p, 


tu 





H 3 

to <U 


tu (U 


1 — 1 

II 

11 

II 

II 

O "O 

3 D, 

rH 

S (X 


Xi 





P p: 

O 6 

d 

B 


to 

(H 

iH 

1 — 1 

1— I 

+J 0) 

o td 

-o 

CO td 


> 

o 

o 

o 

o 

c pH 

to Q 

c 

E-i Q 


o 





o 

•H 

tu 

O 


s 





o 

> 

Ph 

S 


This data was obtained from Reference 12 
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CHAPTER V 
CONCLUSIONS 

A movable mass control system has been conceived to 
convert the tumbling motions of a disabled vehicle into 
simple, spin. It has been shown that a control law re- 
lating mass motion to vehicle motion may be formulated 
using Liapunov stability theory. This technique is use- 
ful for designing control systems where the governing 
equations of motion cannot be linearized. 

For a large space station which is tumbling as a 
result of a collision with a shuttle orbiter , it has been, 
shown that the movable mass system is capable of de- 
tumbling the space station within a period of two hours 
for the assumptions used. This is accomplished using 
a control mass weight of one percent of the vehicle 
weight and a mass amplitude of approximately three meters. 

The following conclusions were drawn concerning 
control system design: 

1. The mass track should be placed as far as 

possible from the vehicle center of mass and 
oriented parallel to the maximum moment of 
inertia axis. 

The control mass weight should be selected 
as large as possible, consistent with peak 
force, power, and energy limitations. 


2 . 


3. The performance of the control system may be 


improved by allowing the mass to travel with a 
larger amplitude. 

From the various cases considered, it can be concluded 
that the parameter selection procedure employed is valid 
and provides realistic estimates of parameter values. 

As an additional application, it has been shown 
that the movable mass control system may be used as a 
wobble damper for a vehicle with an artificial-g mode. 

It has been determined that the performance for this 
system is better than competitive passive devices for 
this application. 

In summary, the movable mass control system can 
effectively damp out unwanted tumbling or nutation 
motions of a large space vehicle with modest sensor, 
force, power,, and energy requirements. 
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